Possible Constraints on the Time Variation of the Fine Structure Constant from 

Cosmic Microwave Background Data 
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The formation of the cosmic microwave background radiation (CMBR) provides a very powerful 
probe of the early universe at the epoch of recombination. Specifically, it is possible to constrain the 
variation of fundamental physical constants in the early universe. We have calculated the effect of a 
varying electromagnetic coupling constant (a) on the CMBR and find that new satellite experiments 
should provide a tight constraint on the value of a at recombination which is complementary to 
existing constraints. An estimate of the obtainable precision is \a/a\ < 7 x 10~ 13 y _1 in a realistic 
experiment. 
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I. INTRODUCTION 

A fundamental question in physics is whether or not 
the physical constants are actually constant. Some uni- 
fying physical theories such as superstring theories do in 
fact suggest that the physical "fine structure" constants 
change in time Q. It is therefore of considerable im- 
portance to find methods of detecting a possible time 
evolution these quantities. In the present paper we wish 
to specifically discuss the electromagnetic fine structure 
constant a. Its present value is known quite precisely to 
be | 



a" 1 = (e 2 /^)- 1 = 1/137.0359895(61). 



(1) 



One option for detecting time variation is of course to 
measure its value in the laboratory and constrain its 
time derivative in this way. However, this has the ma- 
jor drawback that even though quite minute changes are 
detectable the time differences are so small that only a 
moderate sized time derivative is detectable. 

Therefore one often turns to other methods. For in- 
stance it is possible to use astrophysical arguments to 
constrain the evolution of a, the most commonly used 
method being to use differential changes in quasar ab- 
sorption lines. This method offers both a long look back 
time (for z ~ 3 one has £/£o — 1/8, assuming a standard 
flat cold dark matter (CDM) cosmology) and the ability 
to detect rather small changes in a. Another possibility 
is to use Big Bang nucleosynthesis, but this method suf- 
fers from the problem that constraints on a are based on 
a specific assumption on how the neutron to proton mass 
difference depends on a. 

As a possible probe that is complementary to all the 
others discussed we investigate the sensitivity of the 
CMBR to changes in a. It is well known that the fluc- 
tuation spectrum of the CMBR is extremely sensitive to 
the physical conditions at recombination |3|] and, using 
inversion technique, it should therefore be possible to de- 
termine the physical parameters at recombination given 
sufficiently good observations. 



The fluctuations are usually described in terms of 
spherical harmonics 



T(M) =Y,a lm Y lm (0,ct>), 



(2) 



where the coefficients are related to C; coefficients by 

Q = (\a lm \ 2 ). (3) 

These fluctuations were first detected in 1992 by the 
COBE satellite (§, but only for I < 20. At such low I the 
power spectrum is almost degenerate in the cosmologi- 
cal parameters and no real constraints are obtainable. In 
the next few years, however, the power spectrum will be 
measured out to I ~ 2500 by two new probes, MAP and 
PLANCK §, and using this data should yield precision 
measurements of the physical parameters at recombina- 
tion. It should also be possible to constrain new exotic 
physics such as non-standard neutrinos or, indeed, 
a change in a. 

In the next sections we discuss the physical conse- 
quences of changing a and calculate an estimate of how 
precisely we can hope to measure such a change with the 
CMBR data. To calculate actual CMBR power spectra 
we have used the CMBFAST package developed by Sel- 
jak and Zaldarriagga ||. Finally, as will be discussed 
later there is also a terrestrial method which offers long 
look back times, namely to use the Oklo natural fission 
reactor in Gabon. This method offers the currently most 
stringent limit on time variation in a. 



II. CONSEQUENCES OF CHANGING a 

Since the formation of the CMBR is based entirely on 
electromagnetic processes, changing the strength of these 
interactions is bound to change the CMBR fluctuation 
spectrum. First of all, it changes the Thomson scatter- 
ing cross section for all interacting particles. Second, it 
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also changes the recombination of hydrogen. The sec- 
ond effect is far more subtle than the first since it also 
involves changing all the energy levels of the hydrogen 
(and helium) atom. In the following, we shall neglect 
the impact on helium and only concentrate on hydrogen. 
Notice that there is also a small secondary effect from the 
change in helium abundance from nucleosynthesis which 
we shall also neglect in the present paper. 

Thomson scattering - By far the most efficient equi- 
libration mechanism for thermalising the photon gas in 
the early universe is Thomson scattering on free electrons 
(not protons since the rate for this process is suppressed 
by a factor ml /nil ~3x 10 -7 ). The fundamental cross 
section for this process is given by M 



6ir nil 



(4) 



meaning that it has a a 2 dependence. 

Recombination - The phenomenon of recombination 
is of paramount importance for the formation of the 
CMBR since photon equilibration is mediated by Thom- 
son scattering on free electrons. Prior to recombination 
the photons are tightly coupled to the electron-baryon 
fluid, whereas subsequent to recombination the photons 
are therefore essentially free particles. Thus the epoch of 
CMBR formation is directly linked to the recombination 
epoch. The recombination of hydrogen has been exten- 
sively studied by many authors and we shall follow the 
treatment by Ma and Bertschingc r [HTjl which is based on 
the earlier treatment by Peebles 

Recombination directly to the ground state is strongly 
prohibited in the early universe since it leads to immedi- 
ate re-ionisation. Instead it proceeds via 2-photon emis- 
sion from the 2s level or via the redshift of Ly-a photons 
out of the line center [l]J| , Putting together these two 
effects with the appropriate recombination coefficients to 
all exited levels, one obtains an equation for the time- 
evolution of the ionisation fraction, x e = n e /nn, with 
respect to conformal time 



dx e 

— — = aC r 
ar 



f3(T b )(l-x e )-n H a i2) (T b )x 



(5) 



where the first term on the right hand side describes col- 
lisional ionisation from the ground state and the second 
describes the recombination rate, a is the cosmologi- 
cal scale factor, normalised so as to be equal to one at 
present, is the total number density of hydrogen nu- 
clei and T b is the baryon temperature. The other factors 
are given in the following way 



a 



(2) 
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B 1 /T b(x (2) 



Bi = 



2tt 

„2 



/2 = 13.6 eV. 



(6) 
(7) 

(8) 
(9) 



The reduction factor C r has been calculated by Peebles 
Jlil and is given by 



A a + A 2s ^i 



A Q +A 2s ^ ls +/3( 2 )(T b 



where 



pW(T b ) =P{T b )e^' T » 



A a = 
X n = 



a 2 A^(l - x e )n H 
317' 



(10) 

(11) 

(12) 
(13) 



All these equations scale quite straightforwardly with a. 
The only thing left is to treat the two-photon process 
where, in the standard case, A2 S -^i s = 8. 22458s -1 JT2| . 
Following Shapiro and Breit [jl3| one finds that this fun- 
damental process has the very steep dependence 



A 



2s^ls oc a 



(14) 



To see how the process of recombination changes with a 
we have plotted the evolution of the ionisation fraction x e 
as a function of redshift for different values of x = a/cto 
in Fig. 1. If a increases interactions become stronger and 
equilibrium is maintained longer, meaning that the final 
ionisation fraction becomes smaller. This is exactly the 
trend seen in Fig. 1. 
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FIG. 1. The ionisation fraction as a function of redshift for 
three different values of x = a/ao- The solid curve is for 
x = 1, the dashed for x — 0.95 and the dotted for x = 1.05. 

The fact that a lot of the parameters entering the re- 
combination equations are extremely sensitive to changes 
in a brings hope that the CMBR spectrum is equally sen- 
sitive to changes in a. This is exactly the case, as will be 
discussed in the next section. In Fig. 2 we have shown the 
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CMBR fluctuation spectrum for a standard CDM model 
with two different values of x. There are seen to be very 
substantial changes, even for a quite small change in x. 




FIG. 2. CMBR fluctuation spectra for two different values 
of x. The spectrum has been normalised to the quadrupole 
fluctuation 6C2. 



Iii = f> + 1) [Q + C i)errOT r 2 (17) 

where C; )er ror represents the experimental error. Follow- 
ing Lopez et al. JTHJ we shall neglect the experimental 
error and only take into account the "error" induced by 
cosmic variance. It can then be shown that the stan- 
dard error in estimating any parameter is of the order 
of ~ (I~ 1 )u- Specifically, if all parameters are allowed 
to vary simultaneously one obtains 

^ * (1-%, (18) 

whereas if all parameters except Oi have been determined, 
it is 

^^(l u y\ (19) 

Using our cosmological model as given above, we have 
calculated the expected precision to which x = a/a 
can be determined. The results of this calculation have 
been shown in Fig. 3 as a function of the maximum mea- 
sured Z- value. We note here that Z max (MAP) ~ 1000 and 
/max (PLANCK) ~ 2500. Fortunately, since the CMBR 
spectrum is very sensitive to changes in x, it seems pos- 
sible to detect changes as small as 10 -3 — 10 -2 even if all 
cosmological parameters must be determined simultane- 
ously. To be on the conservative side we estimate that 
Sx < 10 -2 is a realistic obtainable precision. 



III. CMBR SENSITIVITY TO CHANGES IN a 



The key question is now whether or not it will be pos- 
sible to detect deviations in a relative to the standard 
value. In order to estimate the sensitivity of the CMBR 
data, we use a standard technique for this purpose. Since 
we have no usable data at present we can only provide 
what is called error forecasting [B| . To do this we choose 
an underlying cosmological model (in our case standard 
CDM) and determine how precisely the cosmological pa- 
rameters can be determined. This method has been de- 
scribed in great detail elsewhere |3|,|l4| and we shall not 
go into details. The cosmological model can be described 
by a vector of parameters and in our calculations we work 
with the following set 




500 



1000 



9 = (fl, fi&, A, h, n, N v , t, a). 



(15) 
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Here r is the optical depth due to possible reionisation 
and n is the spectral index. The standard CDM model 
which we choose as our reference is then given by the 
vector 



©cdm = (1, 0.08, 0,0.5, 1,3,0, ao). 



(16) 



FIG. 3. The expected standard error 8x as a function of 
the maximum measured I in a CMBR measurement of x. 
The solid curve assumes that all other parameters are known 
whereas the dotted line assumes no prior knowledge of any 
parameter. The dashed curve is with fib and h held fixed but 
all other parameters allowed to vary. 



The main point is then to calculate the so-called Fisher 
matrix, which is given by 



Of course in the event that all other parameters can be 
determined by other means it should be possible to detect 
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Sx < 10~ 4 . This is surely not within reach in the foresee- 
able future, but it is still interesting to look at what other 
parameters it is most important to determine in order to 
obtain a better constraint on 5x. It turns out that x is 
most degenerate with J7f, and h, and Fig. 3 we have also 
shown the standard error on x assuming that these two 
parameters are held fixed at their fiducial values ft, = 0.5 
and S1& = 0.08. If these two parameters can be deter- 
mined by other means a factor of 3-5 improvement in the 
precision should be possible. A possible determination of 
J7b should come from Big Bang nucleosynthesis (BBN) ar- 
guments. Especially the new measurements of deuterium 
in quasar absorption systems seem very promising in this 
regard jl6) . As for a measurement of the Hubble param- 
eter perhaps the most promising method is to use what 
is called cosmic complementarity, namely the fact that a 
joint use of CMBR measurements and large scale galaxy 
surveys like the Sloan Digital Sky Survey (SDSS) break 
some of the degeneracy in the CMBR measurements and 
allows for a more precise determination of h jlffl. 

IV. DISCUSSION 

To compare with other constraints we convert the 
above constraint on 5x to a constraint on redshift and 
time evolution of a and obtain 

\oT 1 da/dz\ < 9 x 10~ 5 , (20) 

or 

la^da/dt] < 7 x 10" 13 y _1 . (21) 

This number should be compared with the constraints 
coming from other sources. Especially constraints com- 
ing from the line shift of quasar absorption systems have 
been extremely useful in providing constraints on the 
time evolution of a [^8|-po| . The most recent such mea- 
surement is that of Varshalovich et al. , who obtained 

(Aa/a) *~3 < 1.6 x 10~ 4 . (22) 

This would correspond to 

\a' 1 da/dz\ < 6 x 10~ 5 (23) 

or 

\aT l da/db\ < 1.6 x l(r 14 y -1 - (24) 

It should be noted that there is actually a claim that 
time variation in a has been detected from QSO data 
pl[ . Here, a change of 

Aa/a = -1.5 ±0.3 x 10~ 5 (25) 

has been reported. 

It thus seems that the possible constraints on changes 
in a coming from the CMBR data will be almost as good 



as those from QSO absorption systems if a evolves lin- 
early in time, and potentially better if the evolution is 
non-linear. Moreover it is important to have reliable con- 
straints from different epochs in the evolution of our uni- 
verse. 

We also note that it is possible to constrain the evo- 
lution of a using arguments from Big Bang nucleosyn- 
thesis. However, these are much more model dependent 
than those obtainable from CMBR data. Kolb, Perry 
and Walker p2] | found an upper limit of 

\a~ l da/dt\ < 1.5 x UT^y -1 , (26) 

but, as previously mentioned, this is based on a specific 
assumption of how changes in a affect the neutron to 
proton mass ratio, an assumption which is at best uncer- 
tain. 

Finally there are also quite severe constraints com- 
ing from laboratory experiments and other terrestrial 
sources. Presently, the best laboratory limit is that of 
Prestage, Tjoelker and Maleki |^3j who obtained 

la^da/dt] < 3.7 x l{r 14 y -1 - (27) 

The other very interesting terrestrial constraint comes 
from the Oklo natural fission reactor in Gabon p4] , p5| . 
The most recent discussion is that of Damour and Dyson 
p3| who derived the limit 

|a _1 rfa/dt| < 5 x RT^y -1 - (28) 

In conclusion, we have calculated the expected sensitiv- 
ity of CMBR measurements to changes in the electromag- 
netic fine structure constant a. It was found that CMBR 
should provide a constraint which is on the same order 
of magnitude as other known constraints from cosmol- 
ogy and terrestrial sources. Also, this type of constraint 
is completely independent of all other existing limits, a 
fact which makes it very interesting. It should perhaps 
also be noted here that CMBR data can potentially also 
be used to constrain the time variation of other funda- 
mental constants such as ui e , Gf or G n • 

Note added - After this paper had been submitted an- 
other paper by Kaplinghat, Scherrer and Turner |||] on 
the same subject has appeared. Using the same methods 
they reach conclusions very similar to those presented in 
the present paper. 
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